
PCT/GB 200J / 0 0 5 5 4 5 




INVESTOR m PBOPLB 



^ QflSce I 



PRIORITY 
DOCUMENT 

SUBMITTED OR TRANSMITTED IN 
COMPLIANCE WITH RULE 17.1(a) OR (b) 



I, the undersigned, being an officer duly authorised in accordance with Section 74(1) and (4) 
of Ae Deregulation & Contracting Out Act 1994, to sign and issue certificates onbehalf of the 
Comptroller-General, hereby certify that annexed hereto is a true copy of the documents as 
originally filed in connection with the patent application identified therein. 



The Patent Office 
Concept House 
Cardiff Road 
Newport 
South Wales 




In accordance with the Patents (Companies Re-registration) Rules 1982, if a company named 
in this certificate and any accompanying documents has re-registered under the Companies Act 
1980 with the same name as that with which it was registered unmediately before re- 
registration save for the substitution as, or inclusion as, the last part of the name of the words 
"public limited company" or their equivalents in Welsh, references to the name of the company 
in this certificate and any accompanying documents shall be treated as references to the name 
^ with which it is so re-registered. 



ice with the rules, the words "public limited company" may be replaced by p.l.c, 
or PLC. 



tion under the Companies Act does not constitute a new legal entity but merely 
company to certain additional company law rules. 





Signed 



Dated 19 January 2004 



Art V-va 



10-DEC-2002 13S 11 FROM K R BRYER FIND CO 

Patents Form 1/77 ^ 

Th^KAiiS OFFICE 



TO PATENT OFFICE 



Pabmta Act 1377 



Requeisft for grant of a patent 

S-SSSW'^*"'^*^''*'*' iRKEIVED BY FAX 



1/77 



The Patent onke 
CardxirUcsad 
QweotNPgiRH 



4. 



Your refmnee 



P3055-GB 



Patent plication number 
(Jhe Patent Cffice w^Jittto Ahpi 



0229432.0 



^m'i^ 



^. Full name, address anc! postcsode of the or of 
each applicant ^md&Hine. mmmnesi 



Patents ADP tkumber ffyaa knw 

Jf the a^iicant is a cacporaiB bo^, give the 
country/siatd f>f its mearpmlSion 



Fli^ Soltitions LirniJ^/77oo*"o;oo^^^^ 
St Michad^s Business Centre 
Church Street 
Lyme Regis 
Dorset 
DT7 3DB 



Title of ihcinveniion 



Image Display Systeszn 



S» Name of >our agent {8>9u kmfe onej 

-Address for sefvlot" in the United* Kingdom 
to vhich all corr^jran^encc should b<s sent 
^mdudU^ the pati^e) 



K R Biycf & Co 
7 Gay Street 
Bath 

BAt2m 



Patents ADP inmTaer^^wfeKW /a /d>^'^^ Oo'Z- '^iZ^^^ 



If you are deolariDg priority firom one or moie 
earlier patent appt!cation$, giv^ fho Country 
artel the date of filing of the or of each of diese 
earlieir Applications and (g'you kuAwii) the or 
each applioatlod number 



Pri^Htyspp&ca&Mi number 
(ify<mkn<fw&) 



(defy /moadt /y&wj 



It tM& oi^lkation is divided orottiern^e 
derived Bam an e^irileif UK appUcaiion* 
give the mimbw and lbs jBling date of 
the earlier (^dSoatian 



Nimte of earlier appIicaiicHi 



DafeofftUng 
fjg^y /mpnth /year) 



S« Is a smt^ent of bivenuoivhqi and of right 
10 grant of a patent required in support of 
ibis request? f Answer 'res' (fi 

a) any Appihant named bi parts ffmuan litventar, or 

b) fherefsmfnyensorwhotpntKnamatafan 
applicant or 

c) anyaamedappUeani^at^iprpmmBbod^. 
See note fd}) 



Yes 



Patents Fonn 1/77 



lB-DEC-2002 13S11 



FROM K R BRYER RND CO 



TO PATENT EFFICE 



P-04/34 





Condxuifltioa dieeiB of fonn 




Description 29 



Absltact 



10. If you are also filing way of &fi foDowing 
state how many against cadb. item. 



Priority dociimeAts 



TYansIaiion of priority documoats 

Btsaemsat of inventoiship and right 
to grant cf a patent (Patmar form 7/7?) 

Request fi>r prcUimnHry examination 
and search {Patents Form 9/77) 

Reque&l £br substantiye escamiiiation 
^ateno Fam 10/77) 



After an apptieailcn for a patent kOi he&n ftUsd, tke C^mptrvlier of she Patent Office will consider vfh^th^puhiict;iioh or 
communhadott cftJie invention should be prohibited or resirltted under Section 22 of the Patents Act 1977, You y^ill be informed 
(fit is neeessan^t^profubit or restrict your invenfioft in tfiis^ Furthermore, if you live in the United Kinsdom, Section 23 of 
the Patmts Act 1977 stops you from tmpfyin^for a patent obrOad y/Hhoutjprxt getting written permission from the Paleni Qfflce 
mhss an appUt^on has been filed m least 6 weeks btforekand iti the United Kutsdomjor a patent far the same invention and 
eWker no tiirettdon pmhibWngpuMlcadan or etmummiaaiton has been slven, or otty suoh direotlan lias been ravaked. 



^ If you need he^ tofillinthis form crj^ hove aniy questions^ please contact the Patent Office on 0S4S SOOSOS. 

b) WrlteyoUriUiMerS in OapUal ittters udngblaek ink or you mtty^pe thenu 

c) Jfihere is not enough space for all the relevant details on any pari <f thhform^ please car^fue on a separaHe sheet ^paper 
and write **see conb'nuattOn sheet" in the relevant pari(s). Any contwuadan sheet should bo attached to this form* 

d) !f you have answered Yes' Patents Form 7/77ynilneed to befded 

^ Once you havefllled in the form you tnust remember to sign and date 

J) For tieiails of the fee ttnd v/tiys to pay please eosaact the Patent Office, 



Any oihef dooumeots 
(pleasespe ^ ) 



1/We request the of a patent on the basis of diis ai)]>IiMiion. 




Date 

18 December 2002 



1 2. Name and daytime telephone number of 
pcison to contact in the XJnitBd Klngdnm. 



SHELLEY, Mark 
01225 428877 



Paients Form 1/77 



lB-DEC-2002 13S12 FROM K R BRYER AND CO TO PATENT OFFICE P.06/^ 



In Tunnel Advertising System: Interlacing Image 
Display with Digital Dieplay Panels 



Authors: Richard Wh!te, JeffEvemy 



Overview 

This document describes the In Tuunel Advertising System development to increase 
the di^l^y irame rate by Interlacing images across large digital display panels. The 
display panels include LCD, TFT and other technologies. Each individual frame firom 
a iikn sequence is stretched acro5$ multiple digital display panels that are back lit with 
Xenon strobe lamps that flash simultaneously to display the whole image there by 
making the panel size unimportant. The sequence then moves along flashing at tiie 
train windows to create the animation. Since the Image position relative to the panels 
varies the panel borders are not noticeable in a fiame sequence. 




Flash Separafion 
By generating a frame sequence where images span the 
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History 

On poster-based systems, the fiame rate of a system is determined zs a function of Che 
width of the {m^ge size aad the velocity of the passing train. 















Poster 

Spacing 
^1 


Poster i Poster i 
^ Spacing Spacfng ^! 





FfameRaie = Postens per Second 

Frame Rat© = Train Velocity / Poster Spacing 

LTsing typical poster sLses of A2 or A3 with widths of 667mm 573mni, the 
Ttiinimmn train speeds for safe operation (at 22 finames per second) are 46kph and 
53kph. This presents a limitation for use of the system to track sections where the 
train speeds are always over this and wall result in los$ ofdisplay if trains fall below 
this requirement. 



1 
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Method 

Using digital display pmeil^ to display the rnedia, it is possible to display images ihat 
straddle two cr inore panels; By doing so the limitation of the relationship between 
image size and firarae rate is temoved. A$ the train passes the display, images are 
displayed by simultaneously flashing mnliiple panels across which tlie complete 
picture is stretched. 

Without the technical restrictions of limited frame rates, the. image size is detunttined • 
hy the train window size and distance j&om the tunnel wall. This results in typical 
image sizes of about 1 metre high by 1,8 metres wide. Uring fixed posters this would 
require a miniraum train speed of over 1 40 kph, however using digital displ^ panels 
with this interlacing technique each digital display panel stores several images and 
loads them a$ required. 

To maintain 21 jSnames per second equates to a flash separation of 45ms. 



Rash 

Separation 
Time 



t = 45ms 



^ Window SeparatiOT i 




Distance tnavBlled is 450inni at 10 m/s 



Key: ^.^^ LCD notflaslilng 

LCD flashing 

LCD [mage 

Carriage & Window 

At a train speed of 1 0 m/s (36Iq>h) the flash separation would need to be ^bout 
450mm. Since the image size is about LS metres wide flie digital display panels must 
fl ash 3 or 4 times as the train window / image passes. 
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Train Speed 


S^par^tion 


mis 




mm 


1 


3.6 


45.45 


z 


7.2 


90.91 


3 


10.3 


136.35 


4 


14.4 


1S1,82 


5 


18 


227J27 


6 


21.6 


272.73 


7 


25.2 


31S,18 



Train Speeci " 


Ssparation 


m/s 


kph 


mm 


8 


28.8 


353.64 


9 




40S.09 


10 


35 


454.55 


16 


54 


681.82 


20 


72 


900.09 


30 


108 


1363.64 


40 


144 


1818.18 



Creating an animation for slow track sections, a higher £rame rate sequaice can be 
generated and displayed over a shorter distance. 

In addition to being able to operate in slow tunnel sections $11 locations benefit from a 
smoother image sequence with less nicker. 
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Ught Box /LCD Panel 
Image 



Train Canias© & Window 




Time betwaen ^ ^ — ^ ^ 

consscuttvB' flashes 



Standard Imagd Display without interlacing the image display 



Ught Box / LCD PaneT 
Image 




Train Catiage & Window Cut^i^way ^ 




. t1 t2 . t3 



Time between » 4 > n < — » 

consecutive flashes 



LCD P€inel Image Display with Interiacing the Image display 
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Notes for New patent Application 



Defector Algorithm 



Problem 

Assuming a number of cross track sensors, positioned to detect some feature 
on the train by the breaking or nnaking of a beam- 
In general: 

• The sensors may be optiqel or infra-red, 

• The sensors may use a unidiret^tlonel or br-directional (retroreflective) 
beam. 

• The sensors could use laser beams., . 

• The sensors would be arranged to maximise the clarity of 
measurement of the feature on the train e.g. the gap between the 
carnages. 

The objective Is to process the signals from these sensors In adjacent pairs so 
as to generate the timing for the same feature passing each sensor. The 
problem Is that as the accuracy required is In the order of 1mm the 
measurements are frustrated by the presence of interfering objects such as 
pipes or cables which may be unexpectedly attached to the train. 

In diagrammatic signal form. Figure 1 below presents the Ideal situation: 



J 1 

Figure 1 . Ideal signals from sensors 

Where S1, and S2 represent the signals from each of an adjacent pair of 
sensors with respect to time (the ordinate). The trace SI shows the change In 
signal level as. for example the feature on the train passes causing the status 
of the beam to change for a period W, starting at time T1 . An Identical signal 
appears on trace S2 but at a later time, as the sensor is further along the 
track. In this stndghtfonvard case the timings are given by T1 and T2 
respectively. 



SI 
S2 
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In the more general case, however, ttie effect of Interfsring objects is apparent 
as illustrated in Figure 2 below: 



a be 

SI ru — IS 



S2 



w 



d e f 
Jl I LT 



T1 T2 

Figure 2. Typical signals from sensors 

In tills case the correct timings are again marked T1 end T2 but due to the 
efTecta of Interfering objects there ere other transitions of signals S1 and S2. 
These make the correct transitions less dear and the objective e to provide 
an algorithm to successftjlly extract them. 

Solution 

The solution of this problem requires a decision tree to be traversed as the 
sensor signals are generated. An example Is given overieaf in Figure 4: 

It «n be seen that this algorism has Ivro paths, depending on whether S1 is 
Initially on or off. This has the advantage that features causing occlusion of 
the beam worlc equally with features which clear the beam. For example, for a 
train In which the features are the gaps between carriages (see Figure 3 
below) then detection occurs at the start of the train (beam occluded), and at 
the and of each caniage (clearing of beam) at the points marked 'C. The start 
of subsequent carriages Is ignored as ft occurs during the processing of the 
previous feature (the distance between sensors is less than the length of the 
carriage). By this means, a train of n carnages, will enable n + 1 features to be 
detected at similar intervals. 



□ireciion of travel 



:3 L I I i |— 1 I _ i c 



' I I i i I 

c c c c c c 

Figure 3. Illustration of timing points for six carriage train 
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Rgure 4 Example Sensor Processing Algorfthm 

With reference to Figure 2 the algorithm shown in Figure 4 works as follows: 

As SI Is off at the start the algorithm initialises in State 1 . waiting for Si to 
^^..1^ ^.?J'"*o?*^« algomhm proceeds through State 2 (storing an initial 
value for T1 ) to State 3 and waits for time Ts whilst checking SI is still on. Ts. 
however. Is predetemifned to have a value greater than the width of the pulse 
at a- and ae SI becomes off before Ts has elapsed, the algorithm returns to 
State 0 and restarts. 

At point 'b'. the algorithm proceeds again from State 1 through State 2 (storing 
an updated value for T1) to Stat© 3 and as the length of pulse 'b' Is longer 
than Ts the algorithm then proceeds to State 4 and thence to State 5. now 
ignonng signal SI and monitoring signal S2. 

^If ^"1® ^'9°*''"^n^ proceeds through State 6 (storing an Initial vaiue for 
T2) to State 7 and warts fortlme Ts whilst checking S2 is still on. Ts, however, 
IS precfetemilned to have a value greater than the width of the pulse at 'd' and 
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as S2 becomes off before Ts has elapsed, the algorithm returns to State 5 
and waits for S2 to become on again. 

At point 'e', the aigorithm proceeds again firom State 5 through State 6 (storing 
an updated value for T2) to State 7 and as the length of pulse Is longer 
than Ta the algorithm then proceeds to State 15 and generates a flag to 
Indicate that the values for T1 and 12 are valid. 

The algorithm then restarts and lool^s for the next change In 81. It can be 
seen that; 

• The values T1 and T2 are con-eot, provided that an appropriate value 
for Ts Is chosen. This should he sufficient to Ignore interfering objects 
but not preclude the genuine size of the required feature. 

• The algorithm, under the conditions shown in Figure 3 (assuming SI to 
be on when no train is in the funnel) will first take the route through 
States 8 to 14 (detecting tlie front of the train), and thereafter take the 
route through States 1 to 7 (delecting the back of subsequent 
caniages) and finish at the back of the train. Providing the timing points 
as shown in Figure 3. 



Improvements 

The algorithm may be Improved by any combination of tiie following: 

• A timeout period may be added to States 7 and 14 so that if an error 
has caused the misrepresentation of T1 then the algorithm will ignore 
this and restart. 

• Clearly the effects of Interfering object is asymmetric; this can be 
reflected in the algorithm. For example, pulse 'd' is unlikely to happen 
in reality (unless there is a hole in carriage) and moreover pulse 'e' 
may be short if there Is an interferoig object close to the back of the 
carriage. In this case State 7 should be removed. A similar argument 
could be applied to the removal of State 10. 

« The values of Ts could be difR&rent for States 3, 7, 10 and 14. 

V Additional states could be added to ensure that the detection of Si and 
S2 falls within a window after previous detection events to abrogate the 
effects of serious interfering objects (such as a pantograph halfway 
along a carriage roof). 
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In Tunnel Advertising System 
Calibration Process Using Statistical Profiling 
and IHigh Frequency Drift Containment 



Authars: Matthew Davis, Vicld Parsons^ Jeff 
Evemy^ Jo Cauiirey, 
Vince Johnson, Richard White 



Overview 

This document describes the In Tumel Advertising System calibration process. 
Building Statistical Profiles of the trains passing any given system Dvercomes the real 
world variations in train geometiy. 

The variations are caused by a nunib^ of factors: build variadon^ ride height^ wear 
and tear and non-standard additions to The exterior of the train. There are ftirther errors 
introduced through system installation resulting in differences between the measui'ed 
and actual sensor locations and hdghts; these will also aifect the time^ recorded for 
sensor cuts, and hence the velocities generatcdL 

The differences in geometry from train to train, and against the engineering design 
drawings^ cause the image position in e^ch window to deviate fix^m the expected 
steady mid-position. By analysing the data iSrom many train runs passing the system, it 
is possible to build a profile of train runs then calculate a set of calibraiion values that 
best fit the 'typical* profile. DijQferent types of profile can be recognised in real time as 
the train passes the sysiem and through this calibration process these real-world 
variatiions are eliminated. 

Sines the hiraian eye is very sensitive to sudden movement, it has be«i realised that to 
produce the best perceived display a careful balance has to be stmcJc, minimising the 
overall image drift while eotitafning high Jrequency drift oscillations to within 
acceptable limits. 
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System Anatomy: 

The system measures position, velocity and accelerations &>m passing trains using 
sets of cross^tracfc light-beam sensors. When the front of the tram, or back of each 
carnage, passes a sensor beam the posiUon of the train at that instant Js known. 
Velocities are calculated ftom the time taken for the train to ttavel from one sensor to 

, f?"7r*® ^ described as detectors. Similariy accelerations ate 

calculated ftoni sequential detector velocities. 
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Figure 1 

Measured dimensions of a system that used in the calculation of position, velocity 
and acceleration are SeDsor-Spaclng, System-Offset and Sensot Group Gap The 
accuracy of these measured dim<aasxons critically affect the operation of the system 
and form part of the fimdamental elements lhat must ha calibrated to elimiTiate 
inaccuracies in the building of a system. 
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Figure 2 
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# 

grasses 



As the train progrBifes throuBh a sensor group each event (i.e. the front of the train or 
back of a carnage passing a sensor) is known as a senior cut point. Sensor cut points 
generated flom the second sensorin a group onwards, are known as detector cut 
points^ 



— fTTTff" ff >^ 



Ftve Carriage Trafn Approachtug 3ensor Greup 




IstCutPoirtl 

2nd CutPofnt 
3rd Cut Poj'nl 

4th Cut Point 
SihCut Point 

6th Cut Point 
th Cut Point 
Sth CutPofnt 



Sensors 

Cut Potnte as front of train gind end of each carriage pass sensor 
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•1st Cut Point 
2nd Cut Point 



Sonsor Signal Output as trairi passes 
Figures 



lB-DEC-2^2 13*15 FROM K R BRYER RND 03 



TO POTENT OFFICE 



P. 18/^ 



CalibrtOion Process 

The calibiation process is perfbnned iti two 

• velocity calibration 

• positioiial caKbtatioiL 



Velocity calibration. 

Velocity data from tnany limns h collected. Using this data, the system's sonsor 
(spacing and height) build variations are eliminated. • 

Un0 of best fit 

The -line of best fit' is applied to each set of train data using an appropriate model. 
This may be least squares r^ession, polynomial or other method. An cxainple of this 
IS a Cubic approximation. 

SL'^?r^-1L*?'^ line of best fit the error for each dala point (ie the distance 
trom the data pomts to cubic line) is calculated. 



18 



117 



. ■ * I I 



S 



16 



---44 




I 



J 



i. ... 



©i i 



T 



S 7 e 9 



10 



11 



12 13 M 




Figure 4: Measured velocities for an un-calibrated system, 

Caiibratfng the Sensor Spacing 

The graph of measured velocities for an un-calibrated system, shows an extreme 
exan^pleofa sensor spacing em>r (circle Theciicledpointsbaiongto a 

detector whose measured sensor spacing is smaller than it*s actual spacing This 
produces low velocity readings. 

ll^^^ f "^.^f ^^^^ P^i^t arc tmated as a sensor <.et 

where a single best fit coefficient is calculated and applied to the whole $et. 

Culcuiation af correction coefficient 

Coefficient = m«Aa (i ^ {(model vaXua - mcafiured veloeityj / modeX valua)) 

Where the mean is the error acioss a detector set. 
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Average across marry trains 

The Velocity Calibration is reputed for ramy trains to creal© an average^ correction 
co^cient for each detector. The new spacing valu$$ are moved to the line of best fit 



Slew spacincf » Old a^^acing / (mean (Coef £icient} ) 
This is the mean ooirection coefficient across many trains. 

Individual Cut-Point Calibration 

Having corrected the Sensor Spacing Errors to produce a smooth, -well-fitted data set, 
individual points will still be consistently low or hi^ and off out best fit cwve. By 
analysing many trains these points can be individually calibrated onto the line of best 
fit. The ringed reading at the start and in the middle of the graph are where the front 
of the train enters the fir$t and second group of sensors. 



by 
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Figures 
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These individual cut errors are caused by inaceurate sensor height placing and are 
likely to affect the cut points for front and rear of trains >vifh a sloping cab. 





Pw)5on Error 

Sensor H»}0f^t Hlgber than pravloufi sensor 

Figure 6 

Variatipn in sensor height gives a felse sensor spacing because the front of the train is 
sloping. 

To coxreot these einrors we need to log many trains and look for individual sensor cuts 
that are consistently in exxor. 

Coefficients are applied to each sensor out to bring all values towards the optimal 
•line of best fit\ 

The analysis of the train data i$ performed by computer witli the coefEcients 
programmed into the conltol system and the Brst stage of the calibration is connplete. 
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Image Position Calibration 



Ih this section of calibration the positiott of the image relative to the train window is 
calculated and contiolled to improve the image quality. 

The control system produces the predicted time at which each sensor cut should 
occur. These are compared witii the actual seqasor cut times, which eflfectlvely gives a 
positional error. It is possible to reduce the positiotial error by adjusting the sensor 
groitp positions relative to the system. 

SBits^r str&tch 

Sensor stretch is the co-efficient used for correction of sensor group poallions. This 
will reduce oscilatfons In the image. 

In the first part of the calihration procedure sensor positions are individually 
corrected. The sensor sets are at measured positions in the tunnel, which are relative 
to flie system. As each carriage cut point passes a sensor group the position of the 
image will vary in the window* There is a cumulative error in position that is not seen 
in the previous individual velocity calculations. This produces a difference in the 
predicted image position and the recorded train position. To align the imago in the 
centre of the window, there is a need to adjust sensor group positions, relative to die 
system. 



For given train run the sj^tem records the difierencc in position of image to train. 




Figure? 

Figure 7 shows drift error recorded for a sample train. The zero line is the point where 
the image is central in the window. Points above this line will indicate that the image 
is ahead of the centre of the window, points below the line are bchmd. The red fitted 
linos are appioxknating the gradient of drift readnigs for each sensor set. 

For the example shovwi in Fijjure 7 there is a general movement of the image fiom 
ahead of the central position to behind. This movement is across each sensor group. 
As each feature on the train starts to cross the sensor gioup, the image returns to the 
initial position, ahead of the central point. This lateral movement makes the image 
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'''fif;?/? '5'^ ^^'^"^ ^ quality the difference ftition of ,We to 
minimum). 



In the example shofwn on figure 7, the picture drifte backwards. It is therefore 
nec^sary to ^iaUy increase the distance of the sensors from a prv« *6«yor in the 
group. the first gmup the pivot soosor i& the last sensor position. This is so thai 
when adjustoents are made to the sensor positions, the offset torn the first sroup to 
the start of the system is constant. In all subsequent groups the first sensor is the pivot. 
By keeping the pivot sensor fixed and moving other sensors in the group by a 
correqftonal factor, all drift points win be in error equal to the pivot point In figure 7 

h^.^t^^JS^S'L^J^F^ ^PP"'^^^ -0.05m.Thi5 stretch to the sensor sets 
has the effect of flattening the Imes oti thi$ graph (see figure 9). The average gradient, 
for each sensor set, is taken over all train runs. A diffcxent sensor stretch coeioienf is 
laKcn lor each sensor group. 

Example 1; gradient caiculatioa. 




In an ideal situation the drift should h&. zero. In tkiis example the image drifts lOcm 
over 20^1m of tram distance covered. The distance is in eirxir by lOcm and should 
really mi^asure 20m. In order to bring Sie drift Ihib dovvn to flie horizontal and 
therefore the drift to 2<ao, the gradiait^ taken , 



f * taii-'(0, 1/20-1) = 0.00498 



The mean gradient for all sensor sets over many runs is calculated. This can be 
converted to a stretch by using the equation: 



Stretch ""I - tan(ftv«age gradirait) 

■niia grva a stretch of 0.9950. TMs can be verified by multiplymg the original 

o JSf °" ^^^'^^ ^^'"^S corrected distance of 20m (20. 1 * 0.9950 = 

2to) This works equally fiir negative gradients, ibis gives sensor stretches greater 



End of Example 



iB-DEC-2002 13! 16 FROM K R HWER AND CO 



TO PATENT OFFICE 



P. 23/^ 



After the Stretch 





Sensor gr<f up offsi^t 



FoUowing sensor stretches the drift remains stable. However, as can be seen fix)m 
5aagcpo^tioam.y not be the centre of the WTn^w TlTc^^e 
corrected by shifting the sensor groups by * fixed ofBet amount 



?Sst blTl ^n^^itf^**^.' ^""^ P^^*"« « approximately 

^ : f^^' " Siven an offset of-^0.05m. 

This will move the trigger p^^mt relative the light boxes and so tho picture wiD start 
showmg sooner. The second group is moved relaHve to the first gronp. and so the 
between sensor group gap Is increased. It is likely thai there will have been a 
Tro^^ ^^'^ ^ to accurate to millimetres over ^proximalely 




FIgnre 10 
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Figure 10 shows the drift aftar the sensor group of&ets have been applied. Th^ image 
BOiv centres arornid the zero line at the centre of the tr^in wiruiow. The positional 
error is initially high hut is 5cm or les& for the majority of the nm. 



Filtering and Parameter eon^ot 

t- • * • • 

The' train njms should no-w ^ycsypT^^hti}^ positonal eiror/They can be fuither fine 
tuni-d by adjusting the paratnoteil's ijsed in the application iii the control boxes used 
during an acMjal run. » 

The amomt of error coirected as each new txain.position comes in ia dep^dant on 
several parameters which are sent to the system. These tell the systein the proportions 
in which to use current data vs predicted data» new data vs previous measurements^ 
the point at which to apply any coraections and lai^e error conwtion teims. 

There are a number of methods we combine to find tha parameters that optimise 
performance in a given system. 

a. The first is iteratively, by running the same trains through a simulation 
system and finding v^ues that give the best results- Simulation allows 
identical train runs to be logged with different parabieter seta. The set that 
pioduce:^ the least drift can &en be transfeired to the syatem. 

b. Transforming positiDnal error data from the time domain into tlie lirequsncy 
domain by using a Fourier trajisfonn, and selecting parameter values which 
reduce the hi^ Jiequcncy components to a minimrai. Reduction of high 
frequency noise ensures corrections will be applied gradually producmg a 
smooth picture after qalibration- 

c. By using initial seed values, which have already been deteonined through 
manual calibration. 



Using the filters* Fouriw transfbims and other methods described above it is possible 
to reduce errors such as at the start of the exanaple system. By altering parameters 
used in the real time algorithm train runs with ragidly altering accelerations and other 
behavioural characteristic?^ will have reduced driA on aubsequetil runs as control 
parameters optimise perfomiance for the profile of a given system. 
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Notes for New patent Application 



DifferencB Engine 



Problem 

Assuming updated values for velocity and acceleration there is a naed to 
generate timing pulses to fire tiie Xe lamps in the display boxes In real-time. 

Issues: 

• After Inrttallsation at a fixed position, velocity and acceleration may be 
changed, but the position may not as any discontinuities wif! be \^sible 
and hfghly distracitng. 

» High timing accuracy is required; the display pulses should be accurate 
to within 60JJS or less {corresponding to 1mm at 20ms*^). 

• In order to generate the display pulses, it is necessary to knov^ the 
position of ttie train at ail times. 

• High positional accuracy is required, subsequent frames of each image 
must be positioned to within 1 mm or less to avoid jitter. 

» Timing pulses need to be generated for all the display boxes and for all 
the windows in which the display is to be presented, 

a The positron needs to be available as an Instantaneous value in order 
that extemal processes (for example co-inddence with a sensor beam} 
may use this Information to update acceleration and velocity estimates. 

The objective is to generate the dieplay tfmirtg pulses in real-time from a 
known starting position, assuming updated velocities and acceleration and 
with prior knowledge of the position of the display boxes and the position of 
the windows on the train. 



Solution 

The soluttein of this problem requires Iwo factors: 
1 . The solution of the equatton of motion: 

8 = JJa *dt* + Jv * dt + k (1) 

for position s at a given time t. where a and v are Instantaneous values 
of acceleration (d^s/d^ and velocity (ds/dt) respectively and k is the 
Initial position or ofiset. 
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2. The generation of display pulses using this position Inforniatfon and the 
known posttfon of both wlndov^ and display boxes. 

Equation Qfiyintion 

Assuming the Integration can be perfomned numerically the number of 
calculations can be reduced by rean-anging (1) to give: 

s = + /(v+ /a-dt)»clt + k (2) 

The integrations are definite {from 0 to t) and in numerical terms with discrete 
time units ST may be represented: 

t t 

As 6T remains constant v and a may be appropriately scaled thus: 



(3) 



and: 



\ 



A = a*<n' 



(4) 



(5) 



Assuming that the summation Is perfomied In steps of ST therefore (3) mav 
be simplified to: ' ^ / j 



s= S (V+S A) + k 



(6) 



This calculation may be Implemented using the arohlteoture shown below: 




Figure 1. Implementation of Equation of IMotion 
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The operation can be described as follows: 




marked represent registers of appropriate resolution which 
sntadelay of JT. 



• Acceleration values may be updated at any time. 

• Velocity values may be updated at any Hme. 

• Y®'"® °^ ^ '® speed of the adders and m 
therefore be very small. 

• Resolution of acceleration, velocity and intermedigt© results m 
be arbitrarily defined and Is independent of <5T. 

• Velocity and position values are available to external processi 
at any time. 

Generation of Disoiav Pij| l«^ps 

This information can be used to generate the display pulses by continual 

^mS^"^^.^"" ^"T!"* P®^'^*^" ^9*^^" "^a^ster >K) with that whic 
represents the coincidence of a particular window with a given display box. 

This may be achieved using the algorithm shown In Figure 2 overieaf . 
The nomenclature is as follows: 

• W represents a particular window on the train. 

• B represents a particular display box on the tunnel wail. • 

• tt/(B) represents the next window to pass display box B. 

• Pa represents the position of display box B. 

• PuAB} represents the position of the next window to pass display box B 

• >K represents the instantaneous value of position as calculated above 




Registers gives ttie instantaneous value of position. 
This Implementation has the following advantages: 
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W 4=0 






B = 0 






ctf (B) s 0 




to 


B = B+1 



No 




Is 

B> Number 



YBfi 

■■^^e^ WaftfornewpD$|tton 



B = 0 




Yes 



No 








JiJ B*B + 1 [ 



V) FlashBoxB 






tt/(B) + 1 


45« 1.-.^^ 





No 



^Y^T^ — \ 

B> Number of > 



Ffgure 2. Generation of Display Pulses 
The algonthm operates as follows: 

It is assumed that the algorithm is fnitfalized to State 1 before a train passes 
the system. Thence W end B are Initialized to zero and through States 3, 4. 
and 5 the values of dtf(B) are Initlaiized to zero for afl the display boxes. 
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The algorithm wafts in State 6 until a new position is calculated; it 19 assunri^d 
that this only happens after a train is detected. Then In State 7 the box count 
Is reset to zero and States 8 through 12 are repeated for each box before 
returning to State 5 to wait for the train to move to the next position. 

During State 8 It is determined if the train has moved such that the position of 
the next window coincides with the current display box. If so, a flash is 
generated for that display box and w(B) is set to point to the next window. 

The algorithm ends when the lest window passes the fast box or a time out 
elapses (neither shown In the flowchart), after whfch the algorithm proceeds to 
State 1 and waits for the next train. 

An Important aspect of this algorithm is that it is bounded by the number of 
boxes per control system and is Independent of ffie number of windows, ft is 
only necessary that both the boxes and windows are in the correct sequence. 

More importantfy, because the number of states is bounded, the algorithm 
may be Implemented by a state machine, which operates in finite time; and 
can therefore be Implemented to coincide vwth each step of the position 
generator described previously. Thereby providing the necessary degree of 
accuracy. 



(rtiprovements 

The system as described may be improved by any combination of the 
following: 

* In addition to Instantaneous access, the difference engine (with 
suitable synchronising registers) may be configured to allow 
simultaneous update of any combination of position, velocity and 
aoc^eratlon. This allows any algorithm used to process measurements 
and calculate acceleration and velocity to: 

o Compare positions of events with actual ones; using the error to 
dynamically control the equation of motion (1), 

o Monitor performance and apply meaningful physical limits to 
acceleration and any discontinuities in velocity or position. 

o Avoid disconfinuitfes aitogetfier.by modifying acceleration only; 
using feedbacl^ from both velocity and position errors. 

♦ The design sItowo in Figure 1 has a dear structure whereby the 
acceleration processing is a cascaded element to the velocity 
processing. This may be extended by cascading additionai elements, 
effectively adding third, or higher order integrations to the equation of 
motion (1)- This may be useful for the following reasons: 

o The physical causes (e.g. rate of change of acceleration due to 



application of bral<es etc.) may be known and added to the 
model u^ed to process the measurements, 
o Higher order terms may allow better dynamic control, 
o Higher order terms may assist stabiiify of the overall system . 
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• The algorithm fflustrated In Figure 2 uses values of Pb that clearly 
represent the position of the cflsplay boxes on the tunnel wall. Variation 
of these gives the fbliowing advantages: 

o Allowances may be made for irregular spacing of the display 

boj^, such as to avoid an obstnjctlon on the tunnel wall, 
p The values may be changed to represent the virtual position of 
the box relative to the train. This is Important, for example, when 
part of the display Is situated on a curved section of tunnel. In 
this case equally spaced boxes will appear in different positions, 
relative to their projection onto a curve of larger or smaller 
radius. This &rror can be abrogated by using values of Ps that 
compensate for the difference between arc lengths, 
o By applying a mapping function to the set of values the virtual 
position of any particular display box may be moved. As each 
display conresponds to an individual frame (which is fixed in time 
in the sequence as perceived by the viewer); this pi^vides a 
simple method wrtiereby the position of the display relative to the 
window may be made to follow a predetennlned function of time. 

For example, if we consider the simple case In which we wish to 
move portions of the display (as presented to the viewer) then 
the mapping function fit) is simply the difference between the 
real and virtual posHions and is defined by the following: 

/(t) = p(t)-P (7) 

where t Is the time into the movie, P is the normal position of the 
display with respect to the window, p(t) Is the required position 
as a function oft 

Then, by applying the function to values of Ps to generate new 
values Ps' gives: 

Pb'=/(Tb) + Pb (8) 

where Tb Is the time associated with the frame represented by 
display box B. given by: 

Tb « B / F (8) 
Where F fe the frame rate. 
Combining (8) & (9); 

Pb' = /(B/F) + Pb (10) 

Hence by simply aelecfing an appropriate mapping function /(t). 
and using the "virtual" values Pb In the equatbn of State 8 wa 
can Implement any rrianipulaUon of the Image which may be 
derived or desired. 
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Notes for New patent Application 

Detector Algorithm 

Problem 

Assuming a number of cross tracK sensors, positioned to detect some feature 
on the train by the breaking or making of a beam. 

In general: 

• The sensors may be optlcal'or infra-red. 

• The sensors may use a unidirectional or bi-dlrectlonal (retroreflectivel 
beam. ' 

• The sensors couid use laser beams. 

• The sensors would be arranged to ma>cjmise the clarity of 
measurement of the feature on the train e.g. the gap between the 
carriages. 

The objective la to process the signals from these sensors in adjacent pairs so 
as to generate the timing for the same feature passing each sensor. The 
problem Is that as the accuracy required is in the order of 1mm the 
measurements are fmstrated by the presence of interfering objects such as 
pipes or cables which may be unoipectedly attached to the train. 

In diagrammatic signal fomi. Figure 1 below presents the Ideal situation: 



SI » X- 

-H 



S2 



W 



T1 T2 

Figure 1. Ideal signals from sensors 

Where Si, and S2 represent the signals from each of an adjacent parr of 
sensors with respect to time (the ordinate). The trace SI shows the change In 
signal level as. for example the feature on the train passes causing the status 
of the beam to change for a period W, starting at time T1, An Identical signal 
appears on trace S2 but at a later time, as the sensor is further along the 
track. In this stralghtfonvard case the timings ar^ given by T1 and T2 
respectively. ^ . r 
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In the more general case, however, the effect of Interfering objects fs apparent 
as illustrated in Figure 2 below: 



S1 
S2 



a be 

ru — IS 



W 



d 

n 



e 
r 



f 

ur 



T1 ' ' ' T2 

Figure 2. Typical signals from sensors 



In this case the correct timings are again marked T1 and T2 but due to the 
effects of Interferfng objects there are other transiUonis of signals S1 and S2. 
These make the correct transitions less clear and the objective Is to provide 
an algorithm to successfully extract them. 



Solution 

The solution of this problem requires a decision tree to be traversed as the 
sensor signals are generated. An example is given overleaf in Figure 4: 

it can be seen that this algorithm has two paths, depending on whether SI Is 
initially on or off. This has the advantage that feature^s causing occlusion of 
the beam work equally with features which clear the beam. For example, for a 
ti^in in vifhich the features are the gaps between carriages (see Figure 3 
below) then detection occurs at the start of the train (beam occluded), and at 
the end of each canlage (clearing of beam) at the points marked *C'. The start 
of subsequent canlages Is ignored as it occurs during the processing of the 
previous feature (the distance between sensors is less than the length of the 
carriage). By this means, a train of n carriages, will enable n + 1 features to be 
detected at similar intervals. 



Dfnecfioncf travel 



i i I I i I 

c c c c c c 

Figure 3. Hlustratfon of timing points for six carriage train 
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Ffgure 4 Example Sensor Processing Algorithm 

With reference to Figure 2 the aigorithm shown In Figure 4 works as follows: 

A& S1 is off at the start the algorithm initialises in State 1, waiting for S1 to 
turn on. At point 'a' the algorithm proceeds through State 2 (storing an Initial 
value for T1 } to State 3 and waits for time Ts whilst checiting S1 Is still on. Ts, 
however, is pnedetennined to have a value greater than the width of the pulse 
at 'a* and as S1 becomes off before Ts has alapsed, the algorithm returns to 
State 0 and restart. 

At point 'b'. the algorithm proceeds again from State 1 through State 2 (storing 
an updated value for T1) ta Slate 3 and as the length of pulse 'b' is longer 
than Ts the algorithm then proceeds to State 4 and thence to State 5. now 
ignoring signal SI and monitoring signal S2. 

At point 'd' the algorithm proceeds through State 6 (storing an initial value for 
T2) to Stat© 7 and waits for time Ts whilst checking S2 Is still on. Ts, however. 
Is predetermined to have a value greater than the width of the pulse at 'd' and 
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as S2 becomes on bafore Ts has elapsed, the algorithm returns to State 5 
and waits for S2 to become on again. 

At point 'e', the algorithm proceeds again from State 5 through State S (storing 
an updated value for T2) to State 7 and as the length of pulse is longer 
than Ts the algorithm then proceeds to State 15 and generates a flag to 
Indicate that th^ values for Tl and T2 are valid. 

The algorithm then restarts and looks for the next change in S1. (t can be 
seen that: 

^ The values T1 and T2 are correct, provided that an appropriate value 
for Ts lis chosen. This should be sufficient to ignore interfering objects 
but not preclude the genuine size of the required feature. 

• The algorithm, under the conditions shown in Figure 3 (assuming Si to 
be on when no train Is In the tunnel) will first take the route through 
States 8 to 14 (detecting the front of the train)» and thereafter take the 
route through States 1 to 7 (detecting the back of subsequent 
carriages) and finish at the back of the train. Providing the timing points 
as shown In Figure 3. 



Improvements 

The algorithm may be improved by any combination of the following: 

• A timeout period may be added to States 7 and 14 so that if an error 
has caused the misrepresentation of Tl then the algorithm will Ignore 
this and restart. 

• Cleariy the effects of interiering objects is asymmetric;' this can be 
reflected in the algorithm. For example, pulse 'd' is unlikely to happen 
In reality (unless there is a hole in the carriage) and moreover pulse 'e' 
may be short if there is an Interfering object close to the back of the 
carriage. In this case State 7 should be removed. A similar argument 
could be applied to the removal of State 10. 

• The values of Ts could be different for States 3, 7, 1 D and 14* 

• Additional states could be added to ensure that the detection of SI and 
S2 falls within a window after previous detection events to abrogate the 
effects of serious Interfering objects (such as a pantograph halfway 
along a carriage roof). 
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